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Abstract —Slotted ALOHA can benefit from physical-layer 
network coding (PNC) by decoding one or multiple linear 
combinations of the packets simultaneously transmitted in a 
timeslot, forming a system of linear equations. Different systems 
of linear equations are recovered in different timeslots. A message 
decoder then recovers the original packets of all the users by 
jointly solving multiple systems of linear equations obtained 
over different timeslots. We propose the batched BP decoding 
algorithm that combines belief propagation (BP) and local 
Gaussian elimination. Compared with pure Gaussian elimination 
decoding, our algorithm reduces the decoding complexity from 
cubic to linear function of the number of users. Compared with 
the ordinary BP decoding algorithm for low-density generator- 
matrix codes, our algorithm has better performance and the same 
order of computational complexity. We analyze the performance 
of the batched BP decoding algorithm by generalizing the tree- 
based approach and provide an approach to optimize the system 
performance. 

I. Introduction 

In a wireless multiple-access network operated with the 
slotted ALOHA access protocol, a number of users transmit 
messages to a sink node through a common wireless medium. 
Time is divided into discrete slots and all transmissions start 
at the beginning of a timeslot. Collisions/interferences occur 
when more than one user transmits in the same timeslot (T|. 
Successive interference cancellation (SIC) was proposed for 
slotted ALOHA to resolve collisions so that signals contained 
in collisions can be leveraged to increase throughput (2j- In 
this approach, a number of timeslots are grouped together as 
a frame. Each user aims to deliver at most one packet per 
frame, but it can transmit copies of the same packet in different 
timeslots of the frame. 

To see the essence, suppose we have two users and the 
first user transmits two copies of packet v± in timeslots 1 
and 2 respectively, and the second user transmits one copy 
of packet V 2 in timeslots 2. Since no collision occurs, V\ 
can be correctly decoded by the sink node in timeslot 1. A 
collision occurs in timeslot 2. In SIC, the sink node can use v\ 
decoded from timeslot 1 to cancel the interference in timeslot 
2. This approach can be applied iteratively to cancel more 
interference, in a manner similar to the belief propagation (BP) 
decoding of LT codes over erasure channels. Slotted ALOHA 
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with SIC has been extensively studied based on the AND-OR- 
tree analysis, and optimal designs have been obtained ED-G). 

Physical-layer network coding (PNC) HD (also known 
as compute-and-forward j9j) is recently applied to wireless 
multiple-access network to improve the throughput uni- 
ini- Such multiple-access schemes, called network-coded 
multiple access (NCMA), employ both PNC and multiuser 
decoders at the physical layer to decode one or multiple linear 
combinations of the packets simultaneously transmitted in a 
timeslot. Specifically, Lu, You and Liew m demonstrated by 
a prototype that a PNC decoder may sometimes successfully 
recover linear combinations of the packets when the traditional 
multiuser decoder (MUD) lfl3l that does not make use of PNC 
fails. In the existing works on PNC (or compute-and-forward), 
the decoding of the XOR of the packets of two users has been 
extensively investigated m , iia (see also the overview Ha). 
The decoding of multiple linear combinations over a larger 
alphabet has been studied in J9), fTT) . fl8l . 

In this paper, we consider slotted ALOHA employing PNC 
(and MUD) at the physical layer, called network-coded slotted 
ALOHA (NCSA). We assume that the physical-layer decoder 
at the sink node can reliably recover one or multiple linear 
combinations of the packets transmitted simultaneously in 
one timeslot. Our work in this paper does not depend on 
a specific PNC scheme. Specifically, we consider a A'-user 
NCSA system, where each user has one input packet to be 
delivered over a frame of timeslots. A packet is the smallest 
transmission unit, which cannot be further separated into 
multiple smaller transmission units. But it is allowed to send 
multiple copies of a packet in different timeslots. The number 
of copies, called the degree, is independently sampled from 
a degree distribution. The linear equations decoded by the 
physical layer in a timeslot form a system of linear equations. 
Different systems of linear equations are recovered in different 
timeslots. To recover the input packets of users, a message 
decoder is then required to jointly solve these systems of 
linear equations obtained over different timeslots. Though 
Gaussian elimination can be applied to solve the input packets, 
the computational complexity is 0(K 3 + K 2 T) finite-field 
operations, where T the number of field symbols in a packet. 
In this paper, we study the design of NCSA employing an 
efficient message decoding algorithm. 

With the possibility of decoding more than one linear com¬ 
bination of packets in a timeslot, the coding problem induced 
by NCSA becomes different from that of slotted ALOHA 
with SIC. We will show by an example that the ordinary 


BP decoding algorithm of LT codes over erasure channels 
is not optimal for NCSA. We instead propose a batched BP 
decoding algorithm for NCSA, where Gaussian elimination 
is applied locally to solve the linear system associated with 
each timeslot, and BP is applied between the linear systems 
obtained over different timeslots. The computational complex¬ 
ity of our algorithm is O(KT) finite-field operations, which 
is of the same order as the ordinary BP decoding algorithm. 
We analyze the asymptotic performance of the batched BP 
decoding algorithm when I\ is large by generalizing the tree- 
based approach in lfl9l . We provide an approach to optimize 
the degree distribution based on our analytical results. 

Though the batched BP decoding is similar to the one 
proposed for NCMA E0l . EH . we cannot apply the analysis 
therein. In NCMA, we assume that the number of users is fixed 
but the number of packets to be delivered by each user tends 
to infinity. In NCSA, each user has only one packet while the 
number of users can be large. 

Similar schemes have been developed for random linear 
network coding over finite fields without explicitly considering 
the physical-layer effect, e.g., BATS codes and chunked codes 
(see E 2 I . E3l and the references therein). Here the technique 
for NCSA is different from BATS (or chunked) codes in 
two aspects. First, in BATS codes the degree distribution of 
batches is the parameter to be optimized, while in NCSA the 
degree distribution of the input packets (variable nodes) is 
the parameter to be optimized. Second, the decoding of BATS 
codes only solves the associated linear system of a batch when 
it is uniquely solvable (and hence recovers all the input packets 
involved in a batch), while the decoding of NCSA processes 
the associated linear system of a batch even when it is not 
uniquely solvable. 

In the remainder of this paper, Section|II]formally introduces 
NCSA and presents our main analytical result ffheorem[2}. An 
outline of the proof of the theorem is given in Section [III] An 
example is provided in Section [IV] to demonstrate the degree 
distribution optimization and the numerical results. 

II. Network-Coded Slotted ALOHA 

In this section, we introduce the model of network-coded 
slotted ALOHA (NCSA), the message decoding algorithm and 
the performance analysis results. 

A. Slotted Transmission 

Fix a base field F g with q elements and an integer m > 0. 
Consider a wireless multiple-access network where K source 
nodes (users) deliver information to a sink node through a 
common wireless channel. Each user has one input packet for 
transmission, formulated as a column vector of T symbols in 
the extension field F g m. 

All the users are synchronized to a frame consisting of n 
timeslots of the same duration. The transmission of a packet 
starts at the beginning of a timeslot, and the timeslots are 
long enough for completing the transmission of a packet. Each 
user transmits a number of copies of its input packet within 
the frame. The number of copies transmitted by a user, called 


the degree of the packet, is picked independently according 
to a degree distribution A = (Ai,..., A jf), where D is the 
maximum degree. That is, with probability A^, a user transmits 
d copies of its input packet in d different timeslots chosen 
uniformly at random in the frame. Let A = 1 A (a;) = 

JT AiX 1 and A'(x) = JT iA^x 1-1 . We also call A(x) a degree 
distribution. 

Denote by Vi the input packet of the /-th user. Fix a timeslot. 
Let 0 be the set of indices of the users who transmit a 
packet in this timeslot. The elements in 0 are ordered by 
the natural order of integers. We assume that certain PNC 
scheme is applied, so that the physical-layer decoder of the 
sink node can decode multiple output packets, each being a 
linear combination of v s , s G 0 with coefficients over the base 
field F g . Suppose that B output packets are decoded ( B may 
vary from timeslot to timeslot). The collection of B linear 
combinations can be expressed as 

[«i,...,« B ] = [r s ,s6 0]H, (1) 

where H is a |©| x B full-column-rank matrix over F g , called 
the transfer matrix , and [u s , s G 0] is the matrix formed 
by juxtaposing the vectors v s , where ty comes before v s " 
whenever s' < s". 

Note that in ([TJ, the algebraic operations are over the field 
F g m. We call the set of packets {ui,...,Ub} decoded in 
a timeslot a batch. The cardinality of 0 (the number of 
users transmitting in a timeslot) is called the degree of the 
batch/timeslot. We call the ratio K/n the design rate of 
NCSA. 

Lemma 1. When K/n R as K —> 00 , the degree of a 
timeslot converges to the Poisson distribution with parameter 
A = RA as K —> 00 . 

Proof: This is a special case of Lemma [6] to be proved 
later in this paper. ■ 

Denote by Bid the collection of all the full-column-rank, 
d-row matrices over F g , where we assume that the empty 
matrix, representing the case that nothing is decoded, is an 
element of Bid- For a timeslot of degree d, we suppose that 
the transfer matrix of the batch is H G Bid with probability 
p(H|d). Further, we consider all the users are symmetric so 
that for any d x d permutation matrix P, 

fl(H|d) = p(PH|d). (2) 

The transfer matrices of all timeslots are independently gen¬ 
erated given the degrees of the timeslots. Examples of the 
distribution g will be given in Section HVl 

We say a rate R is achievable by the NCSA system if for 
any e > 0 and all sufficiently large n, at least n(R — e) input 
packets are decoded correctly from the receptions of the n 
timeslots with probability at least 1 — e. 

B. Belief Propagation Decoding 

For multiple access described above, the goal of the sink 
node is to decode as many input packets as possible during a 
frame. From the output packets of the n timeslots decoded by 


the physical layer, the original input packets can be recovered 
by solving the linear equations (|TJ of all the timeslots jointly. 
Gaussian elimination has a complexity 0(K 3 + I\ 2 T) finite- 
field operations when n = ()( K), which makes the decoding 
less efficient when K is large. 

The output packets of all the timeslots collectively can be 
regarded as a low-density generator matrix (LDGM) code. 
Similar to decoding an LT code, which is also a LDGM code, 
we can apply the (ordinary) BP algorithm to decode the output 
packets. In each step of the BP decoding algorithm, an output 
packet of degree one is found, the corresponding input packet 
is decoded, and the decoded input packet is substituted into 
the other output packets in which it is involved. The decoding 
stops when there are no more output packets of degree one. 
However, as we will show in the next example, the ordinary 
BP decoding cannot decode some types of batches efficiently. 
We can actually do better than the ordinary BP decoding with 
little increase of decoding complexity by exploiting the batch 
structure of the output packets. 

For example, consider a batch of two packets u\ and u 2 
formed by 


[ui u 2 \ = [r'i v 2 v 3 v 4 ] 


1 0 
0 1 
1 1 
1 1 


(3) 


to user s. The value of v s is then u : j and hene recovered. The 
algorithm returns the new recovered input packets by searching 
the columns of H with only one non-zero component. 

For a batch with degree d, the complexity of the above 
decoding is 0{d 3 +d 2 T). Suppose that K/n is a constant and 
the maximum degree D does not change with K. Since the 
degree of a batch converge to the Poisson distribution with 
parameter —A (see Lemma Q}. the average complexity of 
decoding a batch is 0(T) finite-field operations. Hence the 
total decoding complexity is 0(KT ) finite-field operations. 

C. Decoding Performance 

For an integer j, denoted by \j\ the set of integers 
{1,.... j}. When j < 0, [j] = 0. For any H £ Pid, define 
7 (H) as the collection of all subsets V of [d — 1] such that 
in the linear system ©, Vi-i(d) can be uniquely solved when 
the values of v r ,r £ i ~ 1 [L] are known. Taking the transfer 
matrix in (OJ as an example, we have 

7 (H) = {{1,3}, {2,3}, {1,2,3}}. 

For a timeslot of degree one, the transfer matrix H is the 
one-by-one matrix with the unity. Then 7 (H) = {0}. For any 
intiger k > 0 , define 

F fc (*)= ^ fl(H|A: + l) £ 2^(1 — tc) fc_|s| . 

h S6 7 (H) 


Suppose that we use the ordinary BP decoding algorithm, and 
when the BP decoding stops, V\ is recovered by processing 
other batches, but v 2 , v 3 and v,\ are not recovered. However, 
if we allow the decoder to solve the linear system ([3]), we can 
further recover v 2 = u 2 — u-\ + v\. The example shows that 
the BP decoding performance can be improved if the linear 
system associated with a timeslot can be solved locally. 

Motivated by the above example, we propose the batched 
BP decoder for the output packets of the physical layer of 
NCSA. The decoder includes multiple iterations. In the i-th 
iteration of the decoding, i = 1 , 2 ,... all the batches are 
processed individually by the following algorithm: Consider 
a batch given in <jT]). Let S' C 0 be the set of indices r such 
that v r is decoded in the previous iterations. When i = 0, 
S = 0. Let ie : 0 —> {1,..., | ©|} be the one-to-one mapping 
preserving the order on 0 , i.e., ie(si) < ie(s 2 ) if and only if 
Si < S 2 - We also write i(s) when 0 is clear from the context. 
The algorithm first substitutes the values of v r , r £ S into <[TJ 
and obtain 

K • • ■, u B \ - [tv, r e S]U^ = [u s , s e 0\S]H i [°\ s l, (4) 

where H 1 -^ is the submatrix of H formed by the rows indexed 
by i[S]. The algorithm then applies Gaussian (Gauss-Iordan) 
elimination on the above linear system so that H 1 [ 0 \‘ s l is 
transformed into the reduced column echelon form H and (Q} 
becomes 

[ui,...,u B ] = K,s e 0\S]H. (5) 

Suppose that the j-th column of H has only one nonzero 
component (which should be one) at the row corresponding 


In other words, T fe ( x ) is the probability that when k + 1 users 
transmitted in a timeslot, the input packets of the user with 
the largest index can be recovered if each of the other users’ 
packet is known with probability x. 

We assume that the maximum degree I? is a contant that 
does not change with I\. The following theorem, proved in the 
next section, tells us the decoding performance of l iterations 
of the batched BP decoder when K is sufficiently large. We 
apply the convention that 0 ° = 1 . 

Theorem 2. Fix real numbers R > 0, e > 0 and an integer 
l > 0. Consider a multiple-access system described above with 
K users and n = \K/R~\ timeslots. Define 

( \ k p~ x \ 

z i = 1 — a ^1 — ^ — r fe (z;_i)J , 

where Zo = 0 and for 1 < i < l 

( \ k p~ x \ 

z i = 1 - A' ( 1 - fc , r fc (^_i)j /A, 

where A = RA. Then for any sufficiently large K, l iterations 
of the batch BP decoder will recover at least K(zf — e) input 
packets with probability at least 1 — exp(— ce 2 K), where c is 
a number independent of K and n. 

Proof: See Section Hill ■ 

Lemma 3. L/A-r) is an increasing function of x. 

Proof: This lemma can be proved by applying ED 
Lemma 13]. ■ 






D. Degree Distribution Optimization 

Theorem [2] induces a general approach to optimize the 
degree distribution A. Let 

f(x; A) = 1 — A' ^1 — X fci TfcOzOj A- 

We have z.i = A ),i = 1,..., l — 1. Suppose that we 

allow l —> oo. The sequence { z ,} is increasing (implied by 
Lemma [3j and converges to the first value x > 0 such that 
f(x; A) = x. For given value of A, 0 < e < 1 and 0 < rj < 1, 
we can optimize the degree distribution A by solving 

max R 

s.t. f{x-,X) > x(l + e), Va; G (0, rj\, (6) 

^2 *Ai = A/f?, ^2 A* = 1, A, ; > 0. 

i i 

Theorem 4 . Denote by R.( A, e) the optimal value of the above 
optimization. Then the rate 

R*( A, e) = R( A, e) ^1 — A ^1 — £ j 

packet per timeslot is achievable for the batched BP decoding 
algorithm. 

Proof: For any 5 > 0, let R = R( A, e) — V~5. We show 
that for sufficiently large K. there exists a degree distribution 
A such that using n < K/R timeslots, the batch BP decoding 
algorithm can recover at least K(rf — \/~5) input packets with 
high probability, where 

->' = ( 1 - A (i-E^r l( ,))). 

That is the code has a rate at least f?*(A, e) — S packet per 
timeslot. 

Let n = \K/R( A, e)]. For the degree distribution A achiev¬ 
ing R( A, e) in (| 6 }, we know by Theorem [2] that at K(z* — \fd) 
input packets can be recovered with high probability. We know 
that the sequence {zi} converges to a value larger than r/. Then 
there exists a sufficiently large l such that zi _ i > r/. Thus, 
z* > rj*. The proof is completed. ■ 

III. Performance Analysis 

We generalize the tree-based approach fl9l to analyze the 
performance of the batched BP decoder and prove Theorem U 

A. Decoding Graph 

The relation between the input packets and the timeslots can 
be represented by a random Tanner graph G, where the input 
packets are represented by the variable nodes, and timeslots 
are represented by the check nodes. We henceforth equate 
a variable node with the corresponding input packet, and a 
check node with the corresponding timeslot. There exists an 
edge between a variable node and a check node if and only if 
the corresponding input packet is transmitted in the timeslot. 
Associated with each check node is a random transfer matrix 


H. For given degree d of the timeslot, the distribution of H 
is g(-\d). 

The /-neighborhood of a variable node v, denoted by Gi(v), 
is the subgraph of G that includes all the nodes with distance 
less than or equal to / from variable node v , as well as all the 
edges involved. Since Gi(v) has the same distribution for all 
variable node v, we denote by Gi the generic random graph 
with the same distribution as Gi(v). After l iterations of the 
batched BP decoding, whether or not a variable node v is 
decoded is determined by its 2 /-neighborhood. 

Motived by the tree-based approach, in the remainder of 
this section, we first analyze the decodable probability of the 
root node of a random tree, and then show that the decoding 
performance of G 21 is similar to that of the tree. The proof of 
Theorem [2] is then completed by a martingale argument. 

B. Tree Analysis 

Fix two degree distributions a(x) and (3(x). Let Ti be a tree 
of / + 1 levels. The root of the tree is at level 0 and the leaves 
are at level /. Each node at an even level is a variable node, and 
each node at an odd level is a check node. The probability that 
the root node has i children is A,. Except for the root node, 
all the other variable nodes have i children with probability 
a,;. All the check node has i children with probability 3 t . An 
instance of X 4 is shown in Fig. Q] 

Lemma 5. Let 2 ;* be the probability that the root variable 
node is decodable by applying the batched BP decoding on 
T 21 . We have 

x* = 1 - A (1 — ]T fc Pk^k(xi-l)) , 
where xq = 0 and for 1 < i < /, 

Xi = 1 - a (1 - . 

Proof: Denote by y, the probability that a check node 
at level 2(1 — i) + 1 can recover its parent variable node by 
solving the associated linear system of this check node with 
possibly the knowledge of its children variable nodes. We have 
x* = 1—A(1 —yi). Suppose that a variable node at level 2(1—i) 
is decodable by at least one of its children check node with 
probability Xi, 0 < i < l. We have x j = 1 — a(l — yf} for 
0 < i < l and Xq = 0 . 

Fix a check node c at level 2(1 — i) + 1. With probability 
<?(H|fc+l)/3fc, the check node has k children and the associated 
linear system has H as the transfer matrix. We permutate the 
rows of H such that the last row of H corresponds to the 
parent variable node. By Q, the permutation does not change 
the distribution g(H|fc + 1). Index the k children by 1 ,k. 
Using Gaussian elimination in the batched BP decoder, the 
parent variable node of check node c can be recovered if and 
only if for certain S £ 7 (H), all the children variable nodes 
indices by S are decodable. Therefore, the probability that the 
parent variable node of c is decodable is X)se 7 (H) ~ 

Xi -for transfer matrix H. Considering all the possible 
transfer matrices, we have yi = ^k^k(xi-i). The proof is 

completed by Xi = Xi. ■ 






We prove the following stronger result than Lemma Q] 

Lemma 6. Suppose that K/n —y R as K —> oo. Fix a timeslot 
t and an integer k > 0. Under the condition that a fixed set of 
k users do not transmit at timeslot t, the degree of timeslot t 
converges to the Poisson distribution with parameter A = RA 
as K —> oo. 


Proof: Let 0 be the set of users that do not transmit 
at timeslot t. For each user that is not in 0, the probability 
that this user transmits a packet at timeslot t is A jn, when n 
is larger than D. Therefore, the degree of timeslot t follows 
a binomial distribution with parameter (K — k,A/n), which 
converges to the Poisson distribution with parameter RA when 
K > oo. ■ 

For a positive integer L, let ez, = 1 — Sd=o ^~7i \—■ We are 
interested in the following instances of a and /3 


ot{x) = = 


1 


-4 


X k e~ x 
k\ ' 


(7) 


Let QifL) be the set of trees of l + 1 levels where each check 
node has at most L children and each variable node has at 
most D children. 


We first consider the case that l is even. Suppose that G;_i 
has N leaf check nodes, which are at level l—l of G;. Denote 
by ki the number of children variable nodes of the i-th check 
node at level l — 1 in G;. Since G; £ Qi{L), we have ki < L. 
By Lemma [ 6 ] we have 

N \ki _A 

Pr{G i = Gz|G i _ 1 = G i _ 1 }^[|—, K ^ oo. 

, r£i’ 

1=1 

On the other hand, we have 


Pi{T l = G l \T l - 1 = Gi- 1 } 


( 1 -cl) 


N 


N 

n 

2=1 


X ki e~ x 

kf- 


Therefore, for sufficiently large K, 


Pr{G ; = G,|G ; _! = G ; _!} - Pr{T; = Gz|T;_! = G,_i} 

> (1 — £l) N — 1 — Cl 

> ~(N + l)e L . 

Note that N = 0(LL z / 2 J). 

We then consider the case that l is odd. Suppose that G/_i 
has N leaf variable nodes, which are at level l — l of G/. 
Denote by ki the number of children check nodes of the i-th 
variable node at level l — 1 of G;. We know that ki < D — 1. 
We then have 


Pr{G; = G;|G;_! = Gz_r} 

TT (ki + l)Afc < + l 
E d dA d 

= Pv{T l = Gi\T l - 1 = Gi- 1 }. 

Therefore, for sufficiently large K , Pr{G/ = G;|G;_j = 

C. Proof of Theorem |2] 


Lemma 7. When I\ is sufficiently large, for any G; £ Qi{L), 

Pr{Gz = G/} > Pr{Tz = G;} - ci, L e L , 

where ci^l = 0 (LLV 2 J) anc i ff le degree distributions ofT\ are 
given in 0. 

Proof: We show by induction that 


Now we are ready to prove Theorem [2] We say G; or Ti is 
decodable if its root is decodable by the batched BP decoding 
algorithm. Fix a sufficiently large L. We have 

Pr{G 2 z £ C? 2 i(£) and is decodable} 

> Pr{G is decodable} Pr{G 2 z = G} 

G602i (£) 


Pr{Gz = G,} > Pr{T; = Gz} - cz, L e L , ( 8 ) 

where c^l = 0 (LLV 2 J)_ 

We prove the lemma by induction. When l = 1, Gi and Ti 
follow the same distribution. For l > 1, we have 

Pr{Gz = Gz} = Pr{Gz = G,| Gz_! = Gz-^PrjGz-! = G,_! 

where G;_i is the subgraph of Gz obtained by removing the 
leaf nodes. We assume that 

Pr{Gz_i = Gz„!} > Pr{Tz_! = Gz_ x } - cz_i, L e i; 

for certain function c;_i.£ = 0 (TL( i - 1 )/ 2 J). We then prove 
© with l > 0 for two cases:/ is even and l is odd. 


> Pr{G is decodable} (Pr{T 2 z = G} - 

G eS 2 z(L) 

> Pr{T 2 z is decodable} — e/4 = xf — e/4 > z* 


--- ) 

4| Q2i{L)\’ 

-e/ 2 , 


where the second inequality follows from Lemma [7] and the 
^ast inequality follows that xf —t zf when L oo. 

Let A be the number of variable nodes v with CXiiiv) £ 
Q 21 (L) and decodable. We have E[A] > (zf — e/2)K. For i = 
1,..., K, denote Z. t = G\[vf). Define X, = E[A|Zi,..., Zf[. 
By definition, X, is a Doob’s martingale with Xq = E[A] 
and Xk = A. Since the exposure of a variable node will 
affect the degrees of a constant number of subgraphs G 2 i(v) 
with check node degree < L + 1, we have \X t — | < 













c', a constant does not depend on K. Applying the Azuma- 
Hoeffding Inequality, we have 

Pr{A < E[A] - e/2 K} < exp ■ 

Hence Pr{A > (zf - e)I\} > 1 - exp(-f^). This 
completes the proof of Theorem [2] 


IV. An Example 


In this section, we use an example to illustrate how the 
proposed NCSA scheme works. Here q = 2 and m = 1. Fix 
an integer N > 2. We consider the PNC scheme that has the 
following outputs: i) When one user transmits in a timeslot, 
the packet of the user is decoded; ii) When two to N users 
transmit in a timeslot, one or two binary linear combinations 
of the input packets are decoded; and iii) When more than N 
users transmit in a timeslot, nothing is decoded. 

Taking N = 3 as an example, when one user transmits in 
a timeslot, the transfer matrix is Hi = [1], and g(Hi|l) = 1. 
When two users transmit in a timeslot, the possible transfer 
matrices are 


H, 


,h 22 = 


h 23 = 


Since H 22 and H 23 have the same probability, g(H 2 i|2) + 
2g(H 22 |2) = 1. Now consider that three users transmit in a 
timeslot. Define 



Y 


'i 


'l 

H 31 = 

l 

,h 32 = 

l 

,h 33 = 

1 1 


l 


i 


1 


The possible transfer matrices are given by the row per¬ 
mutations of H 3 i, i = 1,2,3. Note that for two transfer 
matrices that are permutation of each other, they have the same 
probability to occur. Thus we have 

fl(H 31 |3) + 3 ff (H 32 |3) + 6 5 (H 33 |3) = 1. 

We then have 


r 0 (x) = 1, P i. (x) = g(H 21 |2)x + 2 ff (H 22 |2), 
r 2 (x) = 5(H 31 |3)x 2 + s(H 32 |3)(1 + 2x) 
+ff(H 33 |3)(8x — 2x 2 ). 

In general, for a timeslot of d users, we denote by H,n the 
single column transfer matrix of all ones. For transfer matrices 
of two columns, there are three types of rows: [0,1], [1,0] 
and [1,1]. Denote by Hc; 2 (a) a generic transfer matrix with 
a rows of type [0,1] and d — a rows of type [1,0]. Here 0 < 
a < \_d/2\. All the row permutations of Hd 2 (a) are possible 
transfer matrices. 

Denote by Hd 3 (cii,a 2 ) a generic transfer matrix with a\ 
rows of type [0,1], a 2 rows of type [1,0] and d — a\ — a 2 
rows of type [1,1]. Here a 2 > a\ > 0 and a\ + a 2 < d. 


All the row permutations of Hd 3 (ai,a 2 ) are possible transfer 
matrices. Thus, 


Ld/2j 


i = g (n dl \d)+ J2 

a =1 

d —2 d—ai — 1 

+ £ £ 


0,1 = 1 0,2 = 0,1 


d 

ai,a 2 


ff(H d2 (a)|d) 

g{U d3 (a 1 ,a 2 )\d). 


We can then calculate that 


d- 1 

a — 1 


T d -i(x) = ff(H dl |d)x d_1 

Ld/2j 

+ £ g(U d2 (a)\d) 

a =1 
d —2 d—a± — 1 

+ £ £ s(H d 3 (ai,a 2 )|oO 

0,1 = 1 0,2 — 0,1 

d—1 


x a_1 + 


d—1 

a 


d — 1 

a i — 1, a 2 


„d—a —1 


i —02 — 1 


+ 


+ 


dl, £2.2 — 1 

d—1 

01,02 


i—ai — 1 


d — 0 , 1 — 0,2 — 1 ( ry a 


(x ai + x“ 2 - 


,Q,2 „a,i+a2 


Given average degree A of a timeslot, the average number 
of output packets decoded in a timeslot converges to 

U{ A) = E £ rk(H) fl (H|d), 


H £H d 


when K —> oo. The achievable rate of the NCSA system is 
upper bounded by U (A) packets per timeslot. In the case of 
this example, the achievable rate bound is given by 

Ld/2j 

A) = £ 


X d e~ x 


d\ 


d= 1 

d —2 d—a± — 1 

+ 2 £ £ 

ai=l 02—0,1 


g(U dl \d)+2 ^ f )g(U d2 (a)\d) 

a=l ^ ' 

g(U. d 3 (a 1 ,a 2 )\d) . 


d 

oi, a 2 


Note that the upper bound is in general not tight since the 
packets decoded in different timeslots can be the same. 

We solve © for the above example with the results in Fig. [2] 
where we assume the uniform distribution for each possible 
transfer matrices. We also evaluate the corresponding upper 
bound U( A) for comparison. 
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